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Abstract 


We explore the a - effect and the small-scale current helicity, h c = (h 1 ■ rot h'), for the case of 
weakly compressible magnetically driven turbulence that is subjected to the differential rotation. No 
restriction is applied to the amplitude of angular velocity, i.e., the derivations presented are valid for 
an arbitrary Coriolis number, fi* = 2 Qr cor , though the differential rotation itself is assumed to be 
weak. 

The expressions obtained are used to explore the possible distributions of a effect and h c in 
convection zones (CZ) of the solar-type stars. Generally, our theory gives a^ > 0 in the northern 
hemisphere of the Sun and the opposite case in the southern hemisphere. In most cases the h c has the 
opposite sign to . However, we show that in the depth of CZ where the influence of rotation upon 
turbulence (associated with If*) and the radial shear of angular velocity are strong, the distribution 
of a ^ might be drastically different from a classical cos# - dependence, where 8 is colatitude. It 
is shown that cr ^ has a negative sign at the bottom and below of CZ at mid latitudes. There, the 
distribution of h c is also different from cos#, but it does not change its sign with the depth. 

Further, we briefly consider these quantities in the disk geometry. The application of the devel¬ 
oped theory to dynamos in the accretion disk is more restrictive because they usually have a strong 


differential rotation, 


d log Q/d log r 


> 1 . 
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1 Introduction 

It is generally accepted that the a - effect is one of the most important ingredients of the mean-field 
dynamo (Moffat 1978, Parker 1979, Krause & Radler 1980). Although this effect has been known for 
a long time, there is still some debate about its existence if the nonlinear back reaction of magnetic 
field is taken into account (Moffat 1978, Vainshtein & Cattaneo 1992, Rudiger & Kitchatinov 1993, 
Cattaneo & Hughes 1996, Field et al. 1998, Brandenburg 2001, Ossendrijver et al. 2001). Here we 
shall not discuss this issue, but refer the reader to some recent papers, especially those concerned 
with this question (Field et al. 1998, Brandenburg 2001, Ossendrijver et al. 2001). 

It is well known that the basic dynamo mechanism that is responsible for the generation the 
large-scale solar magnetic field is the combined action of helical turbulent motions, a effect, and 
differential rotation. It is the so-called akl (or, more generally, a 2 Q) dynamo. The question as to 
which extent the differential rotation itself can be responsible for maintaining the a was discussed 
recently by Brandenburg (1999), Rudiger & Pipin (2000) and Rudiger et al.(2001). 

There are several reasons for incorporating the differential rotation in the theory of the a effect: 
l)In the akl - type dynamo, the most important component of the a effect is its azimuthal com¬ 
ponent, a In standard theory, varies as cos#, with # being colatitude (Rudiger & Kitchatinov, 


1 




1993). This simple dependence is probably in contradiction with observations because such an a can 
cause an intense magnetic field at the poles (cf. Brandenburg 1994, Riidiger & Brandenburg 1995). 
This restriction is probably less important for the dynamo operating in the whole CZ. Nevertheless, 
further indirect evidence for a having a maximum magnitude at low latitudes come from examining 
the current helicity observations, see Pevtsov et al (1995) and Kuzanyan et al (2000). Their results 
indicate that current helicity has perhaps a maximum at latitudes near ~ 30°. Both the h c and the 
a effect might originate from a common source. This source is well known. It is either stratified or 
compressible turbulence that is subjected to the Coriolis forces associated with the shear flow (either 
rigid or differential rotation). Rigid rotation gives or’’ oc cos 6. The influence of the differential 
rotation upon the turbulent convection could give a more complicated latitudinal dependence of a . 

2) There are some arguments in favor of the change of the sign of a due to the influence of 
differential rotation if it is strong enough (Brandenburg 1999, Rudiger & Pipin 2000). 

Some of the first results on it were presented by Rudiger & Pipin (2000) and by Rudiger et al. 
(2001). The present paper generalizes those results for the fast rotation case and it is important for 
the astrophysical systems where the case fi* > 1 is quite typical. 

All derivations in the paper are made for the case of weakly compressible magnetically driven 
turbulence. The inclusion of the small but finite compressibility (in the sense that the density fluc¬ 
tuations are allowed for) seems necessary for an existence of a in the originally homogeneous and 
isotropic turbulence subjected to the influence of the mean sheared flow. It is known that the density 
and the turbulence intensity stratifications are probably the most important contributions to a in the 
convection zones of late-type stars (Rudiger & Kitchatinov 1993). However, as a first approximations 
we decided to investigate the influence of the differential rotation on the a and h c in ’’simpler” case 
and then to go ahead in case results are promising. 

Here we assume the turbulence to be magnetically driven in the sense that the original (back¬ 
ground) turbulence consists only of magnetic fluctuations and not of velocity fluctuations. This 
is probably a reasonable assumption for accretion discs where turbulence could be induced by the 
magneto-rotational instability (Balbus & Hawley 1991). The situation in convection zones of cool 
stars is different. In solar plasma the energy of the magnetic part of turbulent energy is likely to be in 
equipartition with the hydrodynamic one (Krause & Radler 1980, Vainshtein 1980, Biskamp 1997). 
So, for reliable estimation of a it is necessary to take both the ’’hydrodynamic” and the ’’magnetic” 
parts of turbulence simultaneously into account as was done, for example, in the paper by Field et 
al. (1999). The disadvantage of their computations is that the spectrum of magnetic and hydro- 
dynamic helicity in the background turbulence were prescribed a priori. We leave the derivation of 
a and h c with the ’’magnetic” and ’’hydrodynamic” parts included simultaneously for the stratified 
differentially rotating flows for a future on the reason said in absatz above. 

In this paper we derive both the a and h c . One reason to do this is that the relation between 
these effects has been commonly used as a diagnostic tool for a-effect in solar physics (Seehafer 1990, 
Pevtsov et al 1995, Kuzanyan et al 2000). Another reason to consider both effects at a time is due 
to the fact fact that magnetic part of a is often associated with the small-scale current helicity. The 
opposite sign of a and small-scale current helicity that is claimed by the Keinigs-Seehafer relations 
(Keinigs 1983; Radler & Seehafer 1990) can be understood from magnetic helicity conservation law 
(e.g. Moffat 1978). The density of magnetic helicity has opposite signs on the large and small scales 
as a result of this conservation law. However, the sign of magnetic helicity of the large-scale fields is 
the same as for a effect. 

The paper is organized as follows, section 2 describes the basic equations and approximations 
we use in our derivations. The general expressions for a and h c are given. Section 3 is devoted 
to examining the main effects in different situations. The applications to the Sun are discussed in 
subsection 3.1. The situation in the disk geometry is briefly considered in the subsection 3.2. Finally, 
the last section summarizes and discusses all the findings. 


2 Basic equations 

2.1 Mean-field electrodynamics 

As usual for the mean Held MHD, we assume the approximate scale separation (Moffatt 1978,2000; 
Krause&Radler 1980). Then the generation of the large-scale magnetic field is described by the 
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following approximation of the mean electromotive force (EMF) of fluctuating fields, 


Si — OLijBj + 


(1) 


with 

£ = (u'xh'). (2) 

It is assumed that the large-scale magnetic field is spatially homogeneous. The current hclicity is 
defined as 

h c = (h ■ rot h ) 

All the derivations below are based on the first-order smoothing approximation. This approxima¬ 
tion is justified if either the magnetic Reynolds number, R m , or the Struhal number, St, are much 
less than 1 (Krause & Radler 1980, Moffat 1978). The first condition is not relevant for astrophysical 
systems. Recently, Field et al (1999) reconsidered the applicability of St << 1 to the mean-field 
MHD. They pointed out the experimentally observed fact that St ~ 0.2 — 0.3 in the ordinary hydro¬ 
dynamic turbulence (Pope 1994). In MHD turbulence, the motions could be largely hydrodynamic 
in character for a modest back reaction of magnetic field, Pouquet et al (1976). Although such a St 
is not very small, we may take it to be a small parameter for pertubation procedure. In addition, as 
pointed out by Moffat (2000), St could be expected rather small on the fast rotating astrophysical 
systems because the turbulence there ”is more akin to a field of weakly interacting inertial waves 
whose frequencies are of the order of the angular velocity of the system”. 

The force field maintaining the turbulence should be defined in the comoving frame of reference. 
In this case it does not contain any information about the mean flow and its gradients, and we can 
safely use the original turbulence concept, which is widely accepted in most mean-field derivations 
(Krause & Radler 1980). This turbulence is supposed to exist in the absence of the mean-magnetic 
field and the mean flow. Note, that our derivations is different at this point from computations made 
by, e.g., Blackman(2000), who uses the concept of the original turbulence, however, the equations 
describing the evolution of the fluctuating fields are written in an inertial frame of reference. 

Thus we have to write at first the equations in the new coordinate system. The mean flow is 
defined by 

Ui = WijXj. (3) 

The tilded quantities are defined in the rest coordinate system. Then 

t 2 

Tij(t) = exp Wt = 6ij + Wijt + — WuWij + - (4) 


The comoving coordinates are given by 

Xi = T^Xj. 

The derivatives are transformed after 


_d_ = T -i_d_ 0_ 
dxi lj dxj ’ dt 


dt 


1 VimXn 


d 

dxi 


(5) 

( 6 ) 


so that the velocity field behaves as 


Ui — 1 VimTmlXl T TijUj. 

The induction equation in the rest frame of reference is 

<9B ~ ~ 

—— = rot { u x B — ri rot B } 
dt 1 ' ! 

In the comoving coordinate system it becomes 

dB ■ 

—^ = (uiBj - UjBi)j + r)ABi - r)t(W p i + Wi p )B ifP i 
at 


to the first order in W. 

A continuity equation has the form, 


dp' dp' , 

-Qj- = WimXm-Q— - p divit . 


(7) 

( 8 ) 

(9) 


( 10 ) 
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Starting from the equation of motion 


one finds that 


din 


1 d 


dt 


p dxi 


B 2 


Bj dBi 


a=T P+7T- +—^+^ i + ^g, 


2 po 1 Pop dxj 


( 11 ) 


_ UV2 ) = ~ 2WilUl “ 2tuWlr 


d Uj 
dxndxi 


(■ Su-t(Wii + Wii )) d 
p dxi 

, Bi 8Bi p' 

H-t;-1-9, 

pop OXl p 


p + 


2-po ) 


( 12 ) 


where the acceleration, g, includes contributions due to gravity and the centrifugal force. Dividing 
the magnetic field for the mean and fluctuating parts, 


Bi = Bi + h[ 0) + hi, 


(13) 


where the contribution to the fluctuating magnetic field itself is made by the background magnetic 
fluctuation, h\°\ (’’the original turbulence”) and h[ are magnetic fluctuations caused by the distortion 
of the mean field, Bi. The latter is governed by the following linearized equation: 

= {u'iBj - u'jBi) d + V Ah'i - V t(W p i + Wip)h' itpl . (14) 

In the comoving frame of reference the fluctuating part of the velocity field satisfies the equation: 

f|r - ^V 2 ) u'i = ^Wuu'i - 2tvWi n /I* (15) 

V (jt / OXnC/%1 

(•Sii-tjWu + Wii )) d ( , | | B t dhf | p' 

p dxi po ) pop dxi p 9 ' 


Next, we extract the solid body rotation part from the mean flow applying Wij = Ei p jQ p + Vij, where 
the term Vij is responsible for differential rotation. 

Upon Fourier-transforming and substituting the last expression for shear, W , we write the induc¬ 
tion equation as 

- - dh' 

{—iuj + gk 2 )h'i = i(k ■ B)ui — 2igV p ik p ki— L (16) 

du> 

and the momentum equation as 


(—iui + uk 2 )u'i + 2(eipjClp + Vij)u'j = - ^iki - V(ip)kpJ^ 


Cgcp’ | (B h i0) ) 
P Pop 


T , , , du[ i(k ■ B) ?(oi p 

- 2ivVpikpki— -1- -h] ’ + — gi, 

dui p 0 p p 


(17) 


where we take into account the relation between density and pressure fluctuations, p = C 2 c p' with 
C 2 c being the sound speed in the turbulent medium, V^ p ) = Vi p + V p i. The compressibility effects 
are described by continuity equation, 


- ioop = —ip(k ■ u) - Vimki 


dp' 

dkjji 


(18) 


h (0) 2 

Equations 11611 71181 can be solved using a perturbation procedure for the small parameters max(u 2 , - 

pop 

and shear, Vij. The first parameter controls the compressibility effects. Finite compressibility should 
be taken into account to obtain the non-zero contributions to a effect. The density fluctuations 
described by ITHl allows for both the Archimedian force and magnetic buoyancy. The combined 
action of these forces and the generalized Coriolis forces produce non-zero contributions to the EMF. 
Actually, the final result contains the factor H^ 1 = g/C 2 c . It can be considered a measure of the 
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typical scale-height in the turbulent medium in which the compressibility effects are important. An 
estimation of H c for the solar convection zone gives H c ~ H v , with H p being the pressure scale 
height. This explains the fact that a effect obtained is the same order of magnitude as a from density 
stratification (cf. Rudiger & Kitchatinov 1993). The whole perturbation procedure for the small but 
finite compressibility was described by Kitchatinov & Pipin (1993). 

The spectrum of background magnetic fluctuations is assumed to be stationary and spatially 
homogeneous, 

(h[ 0) (k,ca)hf ) {k' = ^^7r ij( i(k + k')<5(a; + w'), (19) 

where mj = 5ij—kikj/k 2 . Unfortunately the expressions for the a effect and he in terms of integrals of 
spectral functions are very difficult to manage. So, we have to pass to the mixing-length approximation 
(MLT) in final results. To do this, the procedure proposed by Kitchatinov (1990) is used. We put 
(— iui + gk 2 )~ 1 = (—iui + vk 2 ) _1 = T CO r, with r CO r being the typical correlation time of turbulence. 
The magnetic spectrum is approximated by 


B(k,u)~2{h m 2 )6(k-l^.)6{u). 

(20) 

POO poo 

(/i (0)2 )= / / B(k,u;)dkdcj 

Jo Jo 

(21) 


Results obtained for a and h c are given in subsections below. 


2.2 The a effect. 

Even within the MLT approximation a general structure of the a effect resulting from the shear 
contributions is rather complicated, 

a' 3 = e l e 3 s lmn ( < f29gmV n i+f 3 oe k eiV nk gm+(eg)f3ieiV mn ) ( 22 ) 

+2/ 28 (fi ■ g)e l e 3 + 2/i(fi • g)5 13 + e X3l Vi m (fug™ ~ fi 3 e m (e ■ g)) 

+fnE lm3 Vi l gm + 5 l] e lmn {-f 2 g m V n i + f4e k eiV„kg m + (e ■ g)f 3 eiV m n) 

+Jioo e (e • g)e rn V n i + f 9 e e gie m V' + f 2 se d J e e k g n Vmi 

. p l mnj-rri , p 7 cil mnk-t r . Imn ( n i-w rj , n j T ri\ 

+J6e e Vi e m gn + jiee J d e Vkie m g n +£ g m e n (j 22 e V{ + }ie,e J V) j 
+e 3 £ lrnl Vmn(fi9e n gi + fisg n ei - fn(e ■ g)e"ei) + f 2 oe 3 e lm \e ■ g)V im 
+5 3n E lmz Vmn(fi2gi + /n(e • g)ei) + 2/ 2 7 gW + 2f2ag 3 W 

n fig, 7 i\ Imn t t . n i Imi n t t . r i Imi n t t 

-J 26 (g e J + g J e )e eiV mn + J 24 ,e e J g eiV mn + j 2 ig£ e eiV mn 

+E lmi Vi(f sgi +fr(e ■ g)ei) + g j e lmi (f 33 e n e m V ln + f 32 V lm )} 

J Hop^ac 

where all functions with indices are functions of Coriolis number, fi *. Expressions for all of them can 
be found in the Appendix. For the slow rotation case we have 


' 2g i QP ( Cl-g)S i3 g j Q. i 


. 1 ( 2j ilm , ilm x ri ill mx r ) 

+ 5 V3 £ 9 Vlm ~ £ 9lV ™ _ 3 £ 9 Vlm ) 

c-Hm 1 1 2 /t( 0 ) 2 v 

+ i_J» (Vi + f"V„.) - + VUl] 


It should be noted that for this particular case, the solid body rotation part of the a - effect, the 
upper string, can be restored from the terms contributed by shear, the remaining part of the tensor. 
The substitution Vi m = —EimpQ p should be made to accomplish this step. 
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2.3 The current helicity. 

The derivation of the current helicity was made with an additional restriction, for the sake of sim¬ 
plicity, namely, with the assumption that the azimuthal component of the mean field is predominant. 
The final expression for h c is, 


h c = 

i 

+£ 

+tp4£ 

i ii‘k 
+£ 


y3 2 (2ipo(fl ■ g) + e jk (ipigjVki + faeeiVjigk + ips{e • g )eiVjk)) 
B l B k (ip6gjVu + ipj(e ■ s) e iVji + ipsgjVu) - 4>a£ jk gjB l Bteke^Vi 

ljk gjB l BiVki + ip8£ jk £ lrnn B„giBjeiVmk 


£ B n il) 10 g.jBkemVii 


r(h W 


I'cor go pCac 


m 


For the slow rotation case we obtain, 


h. c = 


—e iSk (B i g j B , Vu - 2 g : B l BiVki 


L15 


—B gjVki) + 


2(H-g) n2l r! or {h^ 2 ) 

5 . lcorgOpC ac 


(24) 


(25) 


Again the last term can be restored from the rest by the method mentioned above. Both (24) and 
(26) are in agreement with previous findings reported by Rudiger & Pipin (2000). 


3 Some applications. 

To find the expressions for a and h c for a particular coordinate system the back substitution V t j —> 
Wij — £i p jfl p has to be made, with Wij as a notation for the derivative (the covariant one in the 
general case) of the large-scale flow. Below, we examine only the most important component of the 
a effect, namely, a . 


3.1 The spherical geometry. 

For the case of the spherical geometry we obtain, 


a 

+ ■ 
h c 

+ ■ 


sin 9 8 log 12 
~2 89~ 


„ , , cos 0 sin 2 9 8 log 12 

— cos9fill A - - -- f a i 

2 a log r 

gTcor{h (0) 2 


(fa 2 - fa 1 Sill 2 9 ) 


POpCac 


, , „ * cos 9 sin 2 9 8 log 12 , 

— cos 9ipoQ A --- ». if hi 


sin 9 8 log 12 
~2 89~ 


('•Phi - sin 2 Oiphi) 


8 log r 
gT 2 or(h (0) 2 )B 
I'corpOpC'ac 


2 


(26) 


(27) 


where f al = (/ 3 - U - /e - h + hW, f a2 = {f 3 + h ~ h ~ h ~ f&W, •phi = (•pi A~ip 3 - •P&W, 
'•phi = (tpi + ip 3 + ■06)12*. Its expressions are given in the Appendix. In the slow-rotation limit we 
have 


a 

h c 


212 * 

T5~ 

ir 

— 


' 3cos9 | gin q d log 12 ^ gr e 


.(hW 2 ) 


d6 


pO pC ac 


cos 9 + sin 9 


8 log 12\ gr 2 or {h W 2 )B 2 


89 


l 2 org0pCac 


(28) 

(29) 


For the SCZ 8 log Q./89 < 0.3 and we can conclude that in the slow rotation limit cr’’ > 0 and h c < 0 
for the northern hemisphere, and they have opposite signs for the southern hemisphere. 
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For the fast rotation case we obtain, 


2ft* 


5 8 7r sin 8 


128 


cos 8 sin 9 


d log ft 


+ 


d log ft 
d9 




gr c , 


Slog r 

2 > 


POpCac 


h c ~ 


cos 8 7r sin 8 , 
■ + 


ft* 


Slog ft 

w~ 


64 

^ + 3 sin 2 ( 


3 cos 8 sin 8 


d^cor 


Slog ft 
S log r 
(h m 2 >B 2 


(30) 


(31) 


Icor go pClc 


where we retain the contributions like ~ 1 /ft* in zero order terms relative to the shear because the 
latter is assumed small. For estimation we put 8 = 7 r/ 4 . Near the bottom of the SCZ where we might 
expect ft* 1, the latitudinal shear of the angular velocity is much less than the radial one. For 
the radial shear there we estimate cos 8 sin 2 #S log ft/S log ?’ < 0.3, see Figure 2 below. So for a , we 
can estimate the conditions on the Coriolis number for equipartition between the solid body rotation 
part and terms contributed by shear, 


1 

ft£ 


97 t S log ft 
256 S log r ’ 


ft* ~ 30. 


(32) 


For h c we have, 

flw = \ fta (33) 

These conditions are probably too severe for the SCZ. However the ft„ may well be expected beneath 
CZ, in the overshoot layer. Thus we can conclude that the condition for sign reversal of the h c ( in 
the sense of the theory developed) is hardly fulfilled on the Sun. 


To estimate the and h c in SCZ more accurately we use the model of the solar interior given 
by Stix (1990) and the helioseismology data reported by Kosovichev et al (1997) with an analytical 
fit given by Belvedere et al.(2000). The distribution of h c depends on the spatial distribution of the 
large-scale magnetic field B$ which is known only hypothetically. Then we introduce the so-called 
’’force-free alpha”, ctff = h c /B 2 (cf. Kuzanyan et al. 2000). Next, we use an assumption of the 
equipartition between the energy of the small-scale magnetic field and the energy of the convective 
flows, 


2 ) 

pop 


(34) 


with u c being the rms convective velocity. We use a standard choice of the mixing-length parameter, 
oimlt = 1.6. To explore the situation just beneath the bottom of the convection zone, rb = 0.715 Rq, 
in the overshoot region, we use the analytical fit for the velocity there given by 


u'o = u b (tanh(70(r - n)) + 1), 


(35) 


with u'b being the convective velocity at the bottom of the SCZ. The resulted depth of the overshoot 
layer is ~ 0.02 J?q. The Coriolis number can be estimated by a procedure described by Kitchatinov et 
al.(2000). Inside the overshoot region, the Coriolis number was set equal to that at the bottom of the 
CZ. The overshooting was followed to 0.69 -Rq. Some important quantities for the SCZ , namely, the 
Coriolis number, the rms convective velocity and the shear contributions, cos 8 sin 2 #<91og ft/9 log r 
and sin 9d log ft/<9# are presented in Figure 1. The first inspection of Figure 1 shows that near the 
bottom of SCZ the radial shear makes a positive contribution to the a effect at the near equatorial 
latitude and a negative contribution at mid and high latitudes. The same can be said about h. c . This 
explains the results presented on the Figures 2 and 3. Figure 2 shows the general distributions of 
the and h c in the meridional section(left. panels) and the radial profiles of these quantities at 
different latitudes. 

The latitudinal distributions at different depths are shown by Figure 3. A complicated latitudinal 
dependence of at the bottom of SCZ is evident. Note, that the obtained value of a// is on 
the order of magnitude in good agreement with observation. This parameter was computed from 
a different set of the solar magnetograms by several authors (Pevtsov et al. 1995, Kuzanyan et al. 
2000). All of them estimated a/f at ~ 10 _ 8 [l/m]. 
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Figure 1: The rotation law is shown in the upper left panel. The upper right panel shows the distribution 
of the Coriolis number, Cl* (solid line), and the convective velocity, u c (dashed line), which is measured 
in [cm/s]. The lower panel shows the shear contributions. The term cos 9 sin 2 0<9 log /<9 log r is shown at 
the lower left, and the term sind<91ogfi/<90 is shown at the lower right. The solid lines are correspond to 
the bottom of the SCZ, 0.715-Rq, the large dashed line to 0.847?q, and the small dashed lines to 0.95f?Q. 






0.72 0.80 0.88 0.96 
DISTANCE 



Figure 2: The resulting distribution of a^ is shown in the upper panel and a// is shown in the lower 
panel. The left-hand parts shows the iso-contours in meridional section. The Radial dependencies at 
different latitudes are shown at the right. The solid line corresponds to 15° latitude , the large dashed 
line corresponds to 45° and the small dashed line refers to the pole. The is measured in [cm/s], and 
a// is measured in [1/cm]. 
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(xIO 3 ) a ** (xIO -11 ) 




Figure 3: Latitudinal distributions of a^ (left), and a// (right), respectively. Solid lines show the values 
at the bottom of the SCZ, 0.715-Rq. We choose to multiply a ^ by factor 100 at this level to present 
all the dependencies at a time. Large dashed lines show the dependence at the level of O.84f? 0 , and the 
small dashed lines correspond to the level of 0.95-Rq. As in the previous Figure, the or^ is in [cm/s] and 
otff is in [1/cm]. 


3.2 a** and h c for disks. 

In disks, where the angular velocity is dependent on the radius, we obtain the following expressions 
for the azimuthal component of a and h c , 


h c = g z 


ipoto* - 


fa 2 d log Q. 


n(/l (0) 


2 d log r 1 

ip h 2 d log n 


2 d log r 


gopCl c 
{h (0) 2 )B% 


r 2 
' cor 


Icor MO pCac 


In the slow rotation limit we reproduce the results reported by Rudiger & Pipin (2000), 


(36) 

(37) 


2 g z Q.* 


15 


3 d log 0 

4 d log r 


corjhW 2 ) 
fJ'OpCac 


he — 


z n* 


i + 


d log ft 


d log r J 


r(h W 2 )B\ 

or PQpCac 


In the fast rotation case we have 


/ 1 7-7T dlogfA Tcor{h^ 2 ) 

dz ^20* 128 <91ogr J popCac 

f 1 77T <9 log 12 \ Tcor(h W 2 )S| 

64 d log r J IcorPopCic 


(38) 

(39) 


(40) 

(41) 


In order to explore the relation between a^ and TLcurr with due regard for the different Coriolis 
numbers and the shear <9 log fl/d log r we construct two functions that define the sign of the effects 
in the upper plane of the disk where g z < 0. Namely, the function 

ao = /in* - (42) 
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Figure 4: The contours for functions controlling the signs of a^ and h c , equations M . Negative val¬ 
ues are shown by the dashed lines. The effects under investigation have opposite signs almost everywhere 
except the region S ~ — 1. 


where S = dlogQ/dlogr defines the sign of the a^, and so does the function 

H d = ^ 0 n* (43) 

for h c - Figure 4 shows the iso-contours for and Hd for the different fl* and S. From Figure 
4 we conclude that the signs of the a^ and h c are opposite nearly for all values of the Coriolis 
number and the shear except the region with shear S ~ — 1. Thus, the Keinig-Seehafer’s identity 
works well for the disks. However, it should be pointed out that application of the theory developed 
here is questionable for the fast rotating disks with |5'| > 1. In our derivations we assumed that the 
differential rotation is weak. Though, the applications for the slow rotation results (40,41), and for 
|S| > 1 could be justified in our approximations. 


3.3 The dynamo with the shear-dominated a effect. 

This subsection is devoted to some preliminary results of applications the obtained a - effect to ID aft 
dynamo model. The considered model is similar to those in papers by Kitchatinov et al.(1994) and 
Kitchatinov & Pipin(1998). The detailed derivations can be found there. Here, we only rewrite the 
mathematical formulation of the problem with taking into account new contribution to a effect due 
to the shear. The system of the simplified ID dynamo equations with regards for the non-linearity 
caused by the influence of the large-scale magnetic field on the A effect driving the differential rotation 
is following, 


dn _ Pm 8 
8t sin 3 9 89 


9 / 9 8Q ~ 

sin 2 6 ( sin 3 9-— — flo cos 9H 
' 89 


dA , -o - 

- = ^aB+- 


8 r l 8 


89 L sin 9 89 


sin 9 A 


— -vn(9)^Q— + - 


sin f) 


sin OB 


(44) 


where B is the toroidal magnetic field, A is the potential the poloidal component of the field and 
= 8 log Q/8r is the radial gradient of the angular velocity. The coefficient, 


a = cos 9 \ 1 + C sh (l(9) cos 9 sin 2 8 


}• 


(45) 
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Figure 5: The latitude dependence of a used in the model. The dashed line shows the ’’standard” 
dependence, a oc cos 9. The solid line shows the shear-dominated a, 14511 . with Q* = 20 and Ci = 3 (cf. 
Figure 3). 


where 


C sh = Ci 


U i(n*) 


defines the distribution of the alpha effect. We introduced another coefficient Ci to match the 
distribution of the alpha effect near the bottom of the SCZ found in the subsection 3.1. The system 
(Hill can be considered as a simplified version of the so-called ’’interface dynamo” which might be 
operating at the interface between CZ and the overshoot region (cf Parker 1993, Charbonneau & 
MacGregor 1997). The magnetic field is measured in units of the field corresponding to equipartition 
between the magnetic energy and the mean energy of the convective flows, the shear, Cl, is measured 
in units of fio /-R© , and time - in units of the typical diffusion time, R@/r/x ■ Pm = vt/^t is a magnetic 
Prandtl number and we shall assume P m = 1 in what follows, and 


is the dynamo number. Functions 


Vt 


i p a 


To 


Ti 


15 

32 B 4 

2 

B2 


4B 2 

~ 3(1 + B 2 ) 2 
2 + 3 B 2 _ 1 

2>/(IW 




1 -B 2 
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arctan B 


describe the magnetic feedback reaction upon the a- and A- effects. 

If the influence of the magnetic field is neglected then equations din yields a shear distribution 


Q. = ^ (5 sin 2 6 — 4) . 


(46) 
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One has a negative radial angular velocity gradient at the high latitudes and a positive gradient at 
the equator. This is in agreement with the helioseismology data. 

The distribution of the a -effect in mi is described by a. It depends on the strength of the shear 
and the Coriolis number. The model in subsection 3.1 gives SI* = 20 near the bottom of the SCZ. 
With the shear mi and Ci = 3 we can match the distribution of the a ^ near the bottom of the 
SCZ found in subsection 3.1. Next picture, Figure 5, shows the dependence of a used in simulations. 
The dependence a oc cos # is shown by the dashed line. The solid line shows the a contributed by 
shear. At the middle latitudes the difference between both alphas is evident. 

The system mi) was solved numerically. The initial field was weak (B 1), and did not has any 
particular symmetry with respect to equator. We followed the field dynamics to a steady-state regime 
that does not depend on the initial field. Only such terminal solutions will be considered later. 


The critical dynamo-numbers and some results about the non-linear behavior of mi with a oc 
cos 9 were considered in more detail by Kitchatinov et al.(1994) and Kitchatinov & Pipin(1998). 
Equations mi give oscillatory solutions for negative dynamo numbers. The most interesting solution 
with long-term modulations of magnetic field and the shear was found by Kitchatinov et al.(1994) 
for a slightly supercritical regime, T> = —6000. Here we repeat these simulations for a different type 
of the a - effect. 

Figure 6 shows some results of the numerical simulations. The upper panel demonstrates the 
long-term evolution of the butterfly diagrams of the toroidal field in coordinates time-latitude for 
a oc cos#. This diagram is in agreement with the previous findings in papers by Kitchatinov et 
al.(1994) and Kitchatinov & Pipin(1998). The middle panel on Figure 6 shows the same for an a 
effect dependent on the shear, mi . and with the same V = —6000. There are several differences 
between pictures. The first one is that the second case does not show the long-term oscillations 
despite the non-linear interplay between magnetic field and shear. Next, the butterfly wings are too 
concentrated near the pole. The last difference is that period of the considered dynamo is about 5 
times larger. This effect is likely due to the space separation between a and SI effects on latitude. 
Such an effect was discussed firstly by Deinzer et al. (1974). This latitude separation is also probable 
reason for the polar concentration of magnetic field. The down panel on Figure 6 shows the evolution 
of the toroidal magnetic field with the general type of a, 1461 . and with V = —5500. The period of 
dynamo is 8 times large there than for the ’’usual” a effect and it is comparable with the solar cycle 
period. 

The critical dynamo number for the dynamo with the shear-dominated a effect is about T> = 
—5400. We also tried to explore the limiting case where the a was defined only by the shear con¬ 
tributions, i.e, a has form like a oc cos#sin 2 9SI. Only steady solutions were found for the negative 
dynamo numbers in this case. This issue is also in agreement with analysis in the paper by Deinzer 
et al. (1974), who also found the steady type solution after the space separation between a and Si 
effects approaches to some limit. 
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Figure 6: The butterfly diagrams of the large-scale toroidal magnetic field evolution for different type of 
the a effect distribution. Time was computed by taking r/r = 10 13 cm 2 /s and Rq = 6.96 • 10 10 cm into 
account. 
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4 Conclusions 


Let us summarize our findings. We have explored the a effect and the small-scale current helicity 
for the case of weakly compressible magnetically driven turbulence that is subjected to differential 
rotation. The derivations presented here were made without restriction on the angular velocity 
amplitude, i.e., for an arbitrary Coriolis number. The differential rotation itself assumed to be weak. 
The main finding is that for the fast rotation case the influence of the differential rotation on the a 
effect and h c can be quite strong. Both a and h c become non-monotonic functions of the latitude 
with two maxima in one hemisphere. One maximum is at mid latitude and another is at the pole as 
in standard theory. 

Beneath the solar convection zone, at the northern hemisphere, is likely to have a negative 
sign at mid-latitudes. However, this result is not conclusive. The problem is that for the magnetically 
driven a effect the "standard”, cos0-like contributions quenches for the fast rotation case as fi* -1 , 
and the shear contributions approaches to a constant. Thus the negative sign of the a is the result 
of such a behavior of the competitive contributions. Since we know the hydrodynamic part of the 
a ^ approaching a constant under the fast rotating limit (Rudiger & Kitchatinov 1993), for a more 
reliable estimation of the a effect the hydrodynamic part is should also be incorporated. 

A debate is ongoing in the literature about whether the relation between a and h c , which was 
proposed by Keinig (1983), Radler & Seehafer (1990) and Seehafer (1994), is a fundamental one. It is 
a very important question because this relation can be used in solar physics as a diagnostic tool to for 
defining the a effect parameters in the Solar convection zone. The present computations show that 
the phase relation between the two holds (in the sense Keinigs’ formula) for the slow rotation case 
eas. 1281291381391 . sf. Rudiger & Pipin (2000), and so does for the fast rotation case eqs. I30I31I40I 
mt . The Keinig-Seehafer’s formula is violated for intermediate values of Q* and shear. This is clearly 
seen in Figures 3 and 4. The validity of Keinig-Seehafer’s identity depends largely on the assumed 
stationarity and homogeneity of the turbulent flows. Such an approximation is probably valid for 
the slow rotation case or for the saturated state in the fast rotation limit. In addition, the phase 
relation can be violated due to contributions from kinetic helicity. It is known that the isotropic a 
effect can be expressed by the sum of kinetic and current helicities, e.g., Field et al (1999). Such a 
representation reflects the fact that a can be divided into the ’’hydrodynamic” and ’’magnetic” parts. 
We skip the ’’hydrodynamic” part in the paper. Hence we can not make a conclusive statement about 
the relation a to h c in a general case. 

Aplication of a ^ contributed by shear to simple ID model shows the drastic differences in evo¬ 
lution of the generated large-scale magnetic field compared to the case where a ^ oc cos 0. The shear 
contributions to cause the spatial separation between a and f 1 generation processes. As noted 
by Deinzer et al(1974), such a separation results to increase of the dynamo period. Our simulations 
show that the interface dynamo driven by shear-dominated a effect (and with parameters adjusted 
to match conditions near the bottom or right beneath the SCZ) has a period about 1 diffusion time 
which is comparable with the Solar cycle period. Though, the obtained butterfly diagrams are not 
in agreement with observations. This could be improve by including the contributions from the 
meridional circulation or from the anisotropic turbulent transport effects (Kitchatinov 1990). Unfor¬ 
tunately, at least for the solar case, the solution of the period cycle problem obtained from this ID 
model can not be extended for the distributed dynamo in the spherical shell, because the radial shear 
is weak in the main part of convection zone and the Coriolis number is not strong enough there. This 
conclusion is well illustrated by Figures 2,3. As can be seen, the influence of the differential rotation 
on the latitudinal distribution of a is quite strong but not enough to be the decisive factor. 

Nevertheless, in the spherical shell dynamos the contributions of the differential rotation to the 
a effect can introduce an interesting component in the non-linear feedback reaction of the magnetic 
field on the a effect. The magnetic field is well known to suppress the a. However, the numerical 
simulations show the radial gradient of the angular velocity is also likely to be increased in the 
regions with the strong magnetic fields (cf. Pipin 1999). Then we might well expect the a effect 
contributed by shear will be less suppressed by the Lorentz forces than contributions connected with 
rigid rotation. This conclusion should be checked by computing the shear-dominated alpha effect 
after taking the influence of magnetic field on the turbulence into account. Such an argument also 
should be considered in numerical simulation devoted to the magnetic quenching of alpha effect (cf. 
Ossendrijver et al 2001). 

Another interesting point which could be discussed in the light of the results obtained connects 
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with idea that in fully developed turbulence the hydrodynamic and magnetic parts of a effect cancel 
each other, Vainshtein (1980). For the rigid rotation case, as we noted above, the magnetic part of 
alpha effect is quenched opposite to increase the Coriolis number. However its hydrodynamic part 
tends to the constant (Riidiger & Kitchatinov 1993). Consequently, these counterparts are hardly 
canceled for the fast rotating system. Our paper (as well as paper of Rudiger & Kitchatinov 1993) 
is based on the considering the real sources of the alpha effect in rotating stratified turbulence. 
This differs from publications were the hydrodynamic and magnetic helicity spectrum are defined a 
priori (eg. Vainshtein 1980, Field et al. 1999). In principle, an ideal choice involves computing the 
combination of the hydrodynamic and magnetic parts of the a effect while including compressibility 
effects equally in both the density and intensity stratifications cases and further taking into account 
the differential rotation, and the non-linear magnetic feedback. We hope to take some steps along 
this direction in future. 
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5 Appendix 

The functions used in this paper are as follows, 

fi = -^5 (n* (3 + 2 O* 2 ) — 3 (l + O* 2 ) arctan(0*)) , 

/2 = (O* (l5 + O* 2 ) + ( — 15 — 60* 2 + O* 4 ) arctan(0*)), 

fi = - err - 57 (V(-2205 - 3360O* 2 - 8760* 4 + 3840* 6 + 730* 8 ) 

J 1920* 9 (l + 0* 2 ) V v ’ 

+9(245 + 4550* 2 + 2000 * 4 - 320* 6 - 210 * 8 + O* 10 ) arctan(0*)) , 

f 6 =- -1 -— (0*(—2205 - 33600* 2 - 1011O * 4 + 240O* 6 + 640* 8 ) 

J 1920*®(1 + O* 2 ) ' v ’ 

—9(—245 - 4550* 2 - 2150* 4 + 110 * 6 + 160* 8 ) arctan(O')), 

f 7 = - -1 - 77 ^91350* + 134400* 3 + 33070* 5 - 21480* 7 - 8700* 9 

960O* 9 (l + O* 2 ) V 

-15(609 + 10990* 2 + 4650* 4 - 1210 * 6 - 102O * 8 - 6 O* 10 ) arctan(0*)^ , 

fa = -qt^- 77 fo*(-9135 - 134400* 2 - 20020* 4 + 4068O* 6 

J 960O* 9 (l + O* 2 ) V v 

+14850* 8 ) + 15(609 + 1099O* 2 + 3780* 4 - 2780* 6 
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— 1950* 8 — 290* 10 ) arctan(0*)^ , 

f 9 = -r-t-— ( 475650* + 705600* 3 + 151330*® - 122520* 7 

J 960O* 9 (l + O* 2 ) V 

-3710O* 9 - 15(3171 + 57610* 2 + 22950* 4 - 7510* 6 - 4860* 8 
—300* 10 ) arctan( 0 *)^ , 

/n =- q— - 77 f91350* + 134400* 3 + 21820*® - 31680* 7 - 8850* 9 

960O* 9 (l + 0 * 2 ) V 

+15(—609 - 10990* 2 - 3900* 4 + 2140* 6 + 1190* 8 + 50* 10 ) arctan(0*)) , 

fi2 =-- 57 ( 0*(—9135 - 134400* 2 - 25870* 4 + 26880* 6 

J 960O* 9 (l + O* 2 ) V v 

+810O* 8 ) - 15(—609 - 10990* 2 - 4170* 4 + 1730* 6 + 110O * 8 
+10O* 10 ) arctan(0*)^, 

fi s = (o*(—135 - 840* 2 + 110 * 4 ) + (135 + 1290* 2 + 50* 4 

+ 110 * 6 ) arctan( 0 *)^, 

/l4 = 640* 5 (1 + O* 2 ) (~ 2 ™‘ + 4n ‘ 3 + + (27 + 5f2 * 2 “ 15 ^ 4 

+70* 6 ) arctan(0*)^, 

/is =-n-t- 57 ( 22050* + 33600* 3 + 16860*® + 3720* 7 - 550* 9 

J 1920* 9 (1 + O* 2 ) V 

+9(—245 - 4550* 2 - 290O* 4 - 820* 6 - O * 8 + O* 10 ) arctan(0*)) , 

fn = -—^- 777 — ( 639450* + 1433250* 3 + 809540*® - 178220* 7 

J 960O* 9 (l + O * 2 ) 2 V 

-233310* 9 - 36150* 11 + 15(1 + 0* 2 ) 2 (-4263 - 2450O* 2 
+960O* 4 + 570O* 6 + 150* 8 ) arctan(0*)^ 

fis = - -1 - 77 f91350* + 134400* 3 + 21820*® - 29880* 7 

J 960O* 9 (l + O* 2 ) V 

-8250* 9 + 15(—609 - 1099O* 2 - 390O* 4 + 202O * 6 + 1110 * 8 
+90* 10 ) arctan(0*)^, 

/ 19 =-—i-- 7 - ( 91350* + 196350* 3 + 103470*® - 23690* 7 

J 1920* 9 (1 + O * 2 ) 2 V 

—26340* 9 - 3540* 11 + 15(1 + O* 2 ) 2 (-609 - 2940* 2 + 1250* 4 
+60O* 6 + 20 * 8 ) arctan( 0 *)^ 

/20 =- ^ - 57 ( 150* + 40* 3 - 30*® - 3(5 + 30* 2 - O * 4 

J 320* 5 (1 + O* 2 ) V v 

+ 0 * 6 ) arctan( 0 *)^, 

/22 =- cr- -^ fo*(2205 + 33600* 2 + 16860* 4 + 804O* 6 

1920* 9 (1 + O* 2 ) V 

+2810* 8 ) - 3(735 + 13650* 2 + 870O* 4 + 390O* 6 
+1630* 8 + 130* 10 ) arctan(0*)^ , 
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■( 


2205ft* + 3360ft* 3 + 876ft* 5 - 384ft* 7 - 73ft* 


192ft* 9 (l + ft* 2 ) 

—9(245 + 455ft * 2 + 200ft * 4 - 32ft* 6 - 21ft * 8 + ft* 10 ) arctan(ft*)) , 
/ 2 5 = ft*(3 + 2ft* 2 ) + 3(1 + ft* 2 ) arctan(ft*)^ , 

/ 2 7 = ^ (—3ft* + (3 + ft* 2 ) arctan(ft*)) , 

/ 2 8 = (ft* (l5 + 4ft* 2 ) — 3 (5 + 3ft*") arctan(ft*)) , 


/29 = 


p- ^ft*(75 + 68ft* 2 + ft* 4 ) 


64ft* 5 (l + ft* 2 ) 

-(75 + 93ft* 2 + 17ft* 4 - ft* 6 ) arctan(ft*)) 


l + ft* 2 ) 2 ( 


15435ft* + 35175ft* 3 + 22212ft* 5 - 76ft* 7 


96ft* 9 (l + ft* 2 ) 

-2895ft* 9 - 283ft* 11 - 3(1 + ft* 2 ) 2 (5145 + 3150ft* 2 - 590ft* 4 
—330ft* 6 + 9ft* 8 ) arctan(ft*)^ , 


— ^ft*(105 + 100ft* 2 + 3ft* 4 ) + 3(—35 - 45ft* 


^ 31 64ft* 5 (l + ft* 2 ) 

—9ft* 4 + ft* 6 ) arctan(ft*)^ , 

/ 32 = ^ g (3ft*(—1 + ft* 2 ) + (3 — 2ft* 2 + 3ft* 4 ) arctan(ft*)), 


foci = 


■ ( 10080ft* + 15120ft* 3 + 3431ft* 5 - 2324ft* 7 - 595ft* 


160ft* 8 (l + ft* 2 ) 

+15(—672 - 1232ft* 2 - 505ft* 4 + 137ft* 6 + 85ft* 8 + 3ft* 10 ) arctan(ft*)), 
(ft* (ill+ 25 ft* 2 ) - (ill+ 62 ft* 2 + 7ft* 4 ) arctan(ft’)) , 


foc2 = ~ 


32 ft* 


h, 5 , 10,21 = -/ 2 ,/i 6 = /4,/24 = -|/ 20,/26 = h, fw = -fn- The current helicity was defined with 
the following functions, 

6ft* + 4ft* 3 — 6 (l + ft* 2 ) arctan(ft*) 


4ft 


*5 


-(ft* (33 + 58ft* 2 + 29ft* 4 ))+3(l +ft* 2 ) (ll + ft* 2 ) arctan(ft*) 
48ft* 5 (l + ft* 2 ) 

(ft* (-645 - 532 ft* 2 + 9 ft* 4 ) 


4>o = 
tpi = 

^’ 2 “ 96 ft* 5 (l + ft* 2 ) 

+3 (215 +249ft* 2 + 37ft* 4 +3ft 

ip 3 =-(ft*(75 + 144ft* + 92ft* 2 + 240ft* 3 + 9ft* 4 + 96ft* 5 ) 

961^ f 1 | j \ 

+3(—25 - 48ft* - 39ft* 2 - 96ft* 3 - lift* 4 - 48ft* 5 + 3ft* 6 ) arctan(ft*)) , 
—6ft*+2 (3 + ft* 2 ) arctan(ft*) 


4ft* 


1p5 = 
l />6 = 


2 ft* (15 + ft* 2 - 8 ft* 4 ) + 6 (-5 - 2 ft* 2 +3 ft* 4 ) arctan(ft*) 
12ft* 5 (l + ft* 2 ) ’ 

ft* (27 + 54 ft* 2 + 23 ft* 4 ) -9(l + ft* 2 ) 2 (3 + ft* 2 ) arctan(ft*) 
24ft* 5 (l + ft* 2 ) ’ 
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n* (-15 

1p7 - 


4 n * 2 + 3 fi* 4 ) +3 (5 + 3 2 - IT 4 + O* 6 ) arctan(fi*) 

48ft* 5 (l + ft* 2 ) 


tps 


1p9 


4 > 10 
tphl 


ft* (3 + 2ft* 2 + 3ft* 4 ) +3 (-l + ft* 2 ) (l + ft* 2 ) 2 arctan(ft*) 
24ft* 4 (l + ft* 2 ) 

ft* (15 + 13 ft* 2 ) - 3 (5 + 6 ft* 2 + ft* 4 ) arctan(ft*) 

24 ft* 5 ’ 

ft* (-9 + ft* 2 ) + (9 + 2 ft* 2 + ft* 4 ) arctan(ft*) 

16 ft* 4 ’ 

= -j-i- tt- fft*(—405 — 228ft* 2 + 73ft* 4 ) 

96ft* 4 (l+ ft* 2 ) V v ’ 


+3(135 + 12111* 2 - 11!1* 4 + 3ft* 6 ) arctan(ft*)) , 

(«• (iiT + i S2 n- + 43 n“) 

-3 (39 + 77ft* 2 + 45ft* 4 + 7ft* 6 ) arctan(ft*)) 
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